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Abstract. First principles density functional theory calculations within the GGA+U 
approximation were performed for LiMn 1 _. E M 2 .02, a candidate cathode material for 
lithium-ion batteries, with (x = 0.25, M=Ni, Fe, Co, Mg), to investigate the effect 
of doping on the destabilization of the monoclinic structure relative to the layered 
rhombohcdral structure. A primary motivation of this work was to determine to what 
extent the predictions of the electronically more realistic GGA+U treatment would 
differ from those obtained within the GGA. Several significant qualitative changes 
are found. For the pristine system in the rhombohedral structure, Mn ions show a 
high-spin state within GGA+U, rather than the low spin (metallic) state found in 
GGA. The doped rhombohedral structure is unstable within GGA+U, rather than 
metastable, as in GGA. In the monoclinic structure, the dopant oxidation states are 
the same (trivalent Fe, divalent Co, Ni) in GGA+U and GGA. Co and Ni ions show 
a higher spin state in GGA+U than in GGA. The divalent dopants destabilized the 
monoclinic structure to a greater extent than trivalent dopants within the GGA+U, 
as expected from previous GGA calculations. Overall, our results suggest that the 
relatively close agreement sometimes found between properties calculated within the 
GGA and GGA+U may be misleading, because the underlying electronic behaviors 
may be profoundly different. 



1. Introduction 

Among the most promising cathode materials for secondary lithium-ion batteries [H[2j[3] 
are lithiated transition metal oxides. To model the electronic and structural properties 
of such materials realistically from first principles has been technically challenging, in 
view of the shortcomings of self-consistent-field treatments of the strongly correlated 
electronic states in many transition-metal oxides, at the level of the local-density (or 
local-spin-density) approximation. It has long been known, for example, that local- 
density-functional theory methods typically predict gapless (metallic) electronic spectra 
for insulating oxides pE]. Despite the unrealistic electronic spectra predicted by local 
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density functional theory, integral properties that depend on total energy, which involve 
summations over all occupied electronic states, such as crystal structure parameters, are 
nonetheless often predicted with reasonable accuracy within local spin density functional 
theory. Even the structural properties, however, may not be properly modeled, if the 
local-spin-density approximation treatment does not capture the correct underlying 
electronic structure, e.g., the transition metal oxidation state, at least qualitatively. 

In previous work [5], the GGA was applied to calculate some electronic and 
structural properties of layered LiMn02 with substitutional dopants, an interesting 
class of candidate lithium-ion battery materials. Of particular interest is the ability 
of divalent dopants to promote stabilization of the layered rhombohedral structure by 
depopulating the majority-spin Mn e g states and thereby destabilize the Jahn- Teller- 
distorted monoclinic structure to a greater extent than trivalent dopants. We anticipate 
that a trend from trivalency to divalency occurs from early to late transition element 
solutes. This simple picture was borne out reasonably well (cf. Fig. (9) in our previous 
paper [5]) by GGA calculations. The stabilization of the rhombohedral relative to 
the monoclinic structure showed an increasing trend, although with some scatter, for 
dopants across the first transition series, with the early transition elements adopting a 
more trivalent character, and the late transition elements a more divalent character. 

Although this picture, in which the solute oxidation state is the property of primary 
importance, is appealing in its simplicity, the inability of GGA, in general, to predict 
physically accurate electronic spectra of correlated transition metal oxides raises doubts 
about its accuracy in predictions of spectrally sensitive properties such as oxidation 
states. It therefore seems desirable to see how well the results of the GGA treatment 
hold up in higher level calculations. 

Among higher- level mean-field treatments, the LSDA-SIC [TJ, GW [8] and LSDA+U 
(or GGA+U) [HI [10] methods come into consideration. We have employed in this 
work the GGA+U method, which is a relatively simple extension of the GGA method 
employed in our previous work [5], with essentially the same computational burden. 
LSDA+U improves the electron-energy-loss spectra and parameters characterizing the 
structural stability of NiO [UJ, [9], compared to LSDA. Improved results of LSDA+U 
predictions for some other oxides [121 HH EE] further demonstrate its merit. 

In this work, we have applied the GGA+U to treat several of the dopants of layered 
LiMn02 considered in our earlier work [5J, to determine whether the dopants show 
dramatically different behavior within this treatment. We find electronic structures 
that are considerably different for the undoped rhombohedral structure, and for the Co 
and Ni doped materials with monoclinic structure. In addition to our calculations for 
the monoclinic and rhombohedral structures, results are presented for the (undoped) 
orthorhombic structure, not treated in our previous work. 

The method of computation is briefly reviewed in Sec. II, section III presents 
results, and conclusions are given in Sec. IV. 
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2. Computational details 

2.1. Method 

We employ in this work the projector augmented wave (PAW) method, as implemented 
in the Vienna ab initio simulation package (VASP) [IB] , with the GGA (generalized 
gradient correction) parametrization given by Perdew et al [17] . In the GGA+U 
extension of density functional theory [9l [TO], an effective on-site Coulomb potential for 
the (i-states is added to the GGA [IT]. For comparison, calculations were also done at 
the GGA level, using the PAW implementation. In previous work [5], GGA calculations 
were performed using ultrasoft pseudopotentials. 

The integration over the Brillouin zone is performed using special k-points generated 
with Monkhorst-Pack [1 8 j indices (4,4,4). The optimized geometry is finally calculated 
for higher numbers of k-points using indices (8,8,8) and a high plane-wave energy cutoff 
of 550 eV to ensure good convergence of stress tensor during cell-parameter relaxations. 
For each system both the volume and internal coordinate relaxations are performed 
during the self consistent energy minimization. An improved tetrahedron method has 
been used for the Brillouin-zone integration [19]. We select a Hubbard parameter U of 
6 eV for Mn and for the dopants Fe, and Co and 7 eV for Ni, guided by Zhou etal. [20J. 
The exchange energy J is fixed at 1 eV [21] . 

2.2. Structure 

Calculations are performed for the orthorhombic, monoclinic, and rhombohedral 
structures of LiMn02, whose space groups are Pmmn, C2/m and R3m, respectively. 
(GGA calculations only for the monoclinic and rhombohedral structures were presented 
in our earlier work [5].) The experimental lattice constants for the orthorhombic [22] 
and monoclinic pj structures are taken as initial values, prior to optimization. For the 
hypothetical rhombohedral structure of pure LiMn02, we take as starting parameters 
the lattice constants calculated by Mishra and Ceder [23], and a ferromagnetic 
spin configuration. For the orthorhombic and monoclinic structures, spin-polarized 
calculations have been performed for both ferromagnetic (FM) and antiferromagnetic 
(AFM) arrangements of the Mn spin. The corresponding magnetic structural unit cells 
are constructed as described by Prasad etal. [5]. In the monoclinic structures, the AF3 
ordering proposed by Singh [24"] has been used for the AFM case. Supercells with four 
formula units for the monoclinic, rhombohedral structures and eight formula units for 
orthorhombic structure were employed. 

In calculations for doped materials, we focus on the monoclinic rather than the 
orthorhombic structure, although experimental investigations of the doped orthorhombic 
system have also been done [25] . The reason the monoclinic system was investigated in 
this work is its similarity to the rhombohedral structure, of which it is a distorted version. 
This similarity suggests a way to analyze the relative tendency of different dopants to 
destabilize the monoclinic distortion. With sufficient doping, the orthorhombic structure 
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would also likely transform to rhombohedral [26J. 

3. Results and Discussion 

3.1. Pristine System 

Equilibrium lattice parameters, calculated within the GGA+U approximation, for 
undoped orthorhombic, monoclinic and rhombohedral LiMnO^, are listed in Table I. 
Results obtained in previous work are also given, for comparison. The largest differences 
occur for the rhombohedral structure, compared to previous GGA calculations [23]. The 
overall ground-state is found to be AF orthorhombic, in agreement with earlier LSDA 
calculations [23J. Longer Mn-0 bond lengths are obtained in all cases, compared with 
earlier LSDA or GGA calculations [23], which is a general feature of GGA+U. The 
underestimation of lattice constants (and bond lengths) within the LSDA is well known, 
particularly for metals with high compressibility [17]. The dependence of structural 
parameters on U can be understood from the weakening of covalent bonding with 
increasing correlation, which shifts the distribution of electronic charge density from 
the interstitial to the atomic regions, and results in weak orbital hybridization in the 
presence of strong coulomb correlation [27J. 

The total and partial Mn d-electron density of states (DOS) in the monoclinic 
(AFM) structure is shown in Fig. 1. The curves above and below the energy axis 
represent the majority and minority carrier states, and the vertical dash-dotted line at 
zero on the abscissa indicates the Fermi energy. The top-most panel, which shows the 
total DOS, agrees with our earlier work[28j. Our calculated bandgap, about 1.5 eV, is 
slightly higher than obtained in a previous [29J LSDA+U calculation for this material. 

The lower panels show site-projected densities of states. The Fermi level appears 
between the occupied and unoccupied Jahn- Teller split e g states of Mn. The minority 
states on the Mn site are almost entirely unoccupied. The last two panels show the O-p 
DOS, which appears at the top of the valence band [28]. 

It was noted in earlier work [5l [30] that Mn ions adopt a low spin state in 
rhombohedral LiMn02 within the GGA. In contrast, GGA+U shows a high spin state 
for trivalent Mn in rhombohedral LiMn02, as would be expected from its behavior in 
other materials. Rhombohedral LiMn02, however, is structurally unstable within the 
GGA+U; only by distorting the structure to a lower (e.g., monoclinic) symmetry can it 
be stabilized. By performing calculations as a function of U — J, we find that the low 
spin state, and stabilization of the rhombohedral structure, occurs for U — J less than 
about 2.3 eV. Since undoped rhombohedral LiMn02 is unstable within the GGA+U, 
the analysis performed in our earlier work [5] , in which low-spin rhombohedral LiMn02 
was found stable within the GGA, is not applicable. We therefore confine our doped 
LiMn02 GGA+U calculations in this work to the monoclinic structure. 
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3.2. Doped systems 

Our previous calculations within the GGA suggested that the relative stability of the 
rhombohedral structure may be promoted by the introduction of dopants with trivalent, 
and particularly divalent oxidation states, which suppress the cooperative Jahn- Teller 
effect [El [5] . In this work, we investigated the substituents Ni, Fe, Co and Mg on the Mn 
sublattice in AFM LiMn02, in the presence of coulomb repulsion parameter U. 

Fig. 2 shows the Mn and Ni d-DOS for LiNio.25Mno.75O2. The top-most panel (a) 
shows the majority- and minority-spin DOS for Ni, and panels (b), (c) and (d) show 
the majority and minority DOS for the three Mn ions in the four formula-unit supercell 
used in our calculations. The empty majority e g band in panel (b), in contrast with 
the Mn d-DOS for pure LiMn02 (Fig. 1), reveals the oxidation of the Mn to the 4+ 
state in the presence of the Ni ion. Confirming that the Mn ion represented in panel (b) 
is tetravalent is the surrounding oxygen octahedron, which shows essentially no Jahn- 
Teller distortion, unlike the other Mn ions in the cell. Panels (c) and (d) show the DOS 
of Mn 3+ ions, similar to Mn ions in the undoped system, shown in middle panel of Fig. 
1. The majority Ni d-DOS in panel (a) shows fully occupied t 2g and e g bands while 
minority DOS shows fully occupied t 2g and empty e g bands, consistent with a divalent 
oxidation state. A similar analysis of site-projected DOS, as well as oxygen octahedra, 
for the doped systems LiCoo.25Mno.75O2, LiFeo.25Mno.75O2, and LiMgo.25Mno.75O2 reveals 
that Fe adopts the 3+ oxidation state whereas Co and Mg are 2+. 

In the case of Co, divalent or trivalent oxidation states can be obtained, depending 
on the initial atomic arrangement, prior to relaxation of the internal coordinates. If 
the Jahn- Teller distortion of the oxygen octahedron surrounding one of the Mn ions is 
suppressed in the initial configuration, the relaxed calculation yields a tetravalent state 
for than Mn, and a divalent state of Co; otherwise, all transition ions are trivalent. We 
find that the state in which Co is divalent is lowest in energy, and the state in which Co 
is trivalent is therefore metastable. 



3. 3. Destabilization of monoclinic structure by dopants 

Since the rhombohedral structure of LiMn02 is unstable, within the GGA+U 
approximation, we perform a different analysis from that in our previous work [5] to 
assess the destabilization of the monoclinic structure with doping. We note that the 
layered rhombohedral structure is a special case of the monoclinic structure that occurs 
when 

a/b = -y/3, ccos/? = — a/3, (1) 

where a, 6, and c are the lattice constants and (3 is the monoclinic angle In the 
pristine system, 

h(x = 0) = a/(y/3b) = 1.12, (2) 

When dopants are introduced, the relaxed lattice constants a(x), b(x), and c(x) vary. 
Consequently h decreases, and, at a critical concentration that corresponds to the 
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phase boundary between monoclinic and rhombohedral structures, decreases, perhaps 
discontinuously, to unity. (A two-phase coexistence region is also possible, in principle. 
For the sake of discussion, we also ignore solubility limits on x). 
In Fig. 3, we plot the quantity 

H(x) = [h(x)-h(0)]/[l-h(0)], (3) 

evaluated with lattice constants calculated within the GGA+U approximation. The 
function H{x) is, by construction, zero at x = and unity at the critical concentration, 
x cr . The more effective a dopant is in destabilizing the monoclinic structure relative to 
the rhombohedral structure, the larger the value of H(x). (We note that this analysis 
does not directly yield a prediction of the critical dopant concentration.) As shown in 
the Figure, the divalent dopants Mg, Ni and Co are more effective than the trivalent 
dopant, Fe. This is consistent with trends observed in previous work. 

We have also calculated amplitudes of Jahn- Teller modes for Mn- and dopant 
centered octahedra in the monoclinic structure. Our values are close to the 
corresponding values of our previous work[5j indicating that the dopants reduce the 
Jahn- Teller distortion in the GGA+U also. 

4. Summary and Conclusions 

Little difference was found between GGA+U and GGA calculations of some structural 
properties of undoped and doped LiMn02- Lattice parameters calculated for 
orthorhombic and monoclinic structures of pure LiMn02 within the GGA+U framework 
are in reasonable agreement with earlier work based on the GGA. Consistent with 
previous GGA-based calculations [23J, as well as experiment [31], the orthorhombic 
structure in the AF configuration is found to be the overall ground state in the GGA+U 
approximation. 

Significant qualitative differences between electronic structures calculated in 
GGA+U and GGA were found in some cases, however, which indicate that the 
agreement cited above may be misleading. For rhombohedral LiMn02, GGA+U gives a 
high-spin solution, in contrast to the low spin result found in GGA j5[ [30] . Furthermore, 
the rhombohedral structure in GGA+U is not metastable, as in GGA, but unstable. 
Although it cannot be verified experimentally, because rhombohedral LiMn02 is a 
hypothetical structure, the electronic structure calculated within the GGA for (low- 
) doped or undoped rhombohedral LiMn02 are assumed to be unphysical. To treat 
doping in the present work, we have therefore focused on the monoclinic structure, and 
the changes induced by doping. For Ni and Co dopants, divalent oxidation states are 
found both in GGA and GGA+U, but GGA+U gave a significantly higher state of 
polarization. In the case of Fe dopants, a high spin state is predicted both by GGA and 
GGA+U. 

The greater destabilization of the monoclinic structure induced by divalent than 
by trivalent dopants previously found within the GGA, was found to hold also within 
GGA+U. 
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TABLE I. Lattice parameters of monoclinic (MLA), orthorhombic (ORTHO) and 
rhombohedral (RLA) structures of LiMnC>2 for AF and FM configurations, based 
on GGA+U approximation. The values given in parentheses are results of previous 
LSDA(GGA) calculations [23J and the available experimental work [2j [22] respectively. 



Structure 


a (A) 


b (A) 


c(A) 


P 




internal parameter 


MLA(AF) 


5.49 


2.86 


5.44 


112.87° 




x=0.271, z=0.769 




(5.54) 


(2.77) 


(5.47) 


(116.0°) 




(x=0.271, z=0.762)(Ref. 23) 




(5.44) 


(2.81) 


(5.39) 


(116.0°) 




(x=0.272, z=0.771)(Ref. 2) 


(FM) 


5.48 


2.87 


5.43 


112.85° 




x=0.270, z=0.769 




(5.54) 


(2.82) 


(5.44) 


(116.0°) 




(x=0.270, z=0.763)(Ref. 23) 


ORTHO (AF) 


2.87 


4.60 


5.81 


90.0° 




zr,=0.117, zmt7=0.636 
z =0.139, 0.600 




(2.79) 


(4.69) 


(5.64) 


(90.0°) 


(Zii : 


=0.104, z M n=0.636)(Ref. 20) 
(z =0.132, 0.606) 




(2.81) 


(4.57) 


(5.76) 


(90.0°) 


(,ZLf- 


=0.126, z M „=0.635)(Ref. 22) 


(FM) 


2.87 


4.61 


5.81 


90.0° 




z Li =0.117, z Mn =0.637 
z =0.138, 0.600 




(2.80) 


(4.82) 


(5.60) 


(90.0°) 




=0.116, z Mn =0.638)(Ref. 23) 
(z =0.128, 0.605) 
(z =0.144, 0.602) 


RLA(FM) 


3.02 


3.02 


14.58 


90.0° 




x=0.257 




(2.82) 


(2.82) 


(14.27) 


(90.0)° 




(x=0.255)(Ref. 23) 
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FIGURE CAPTIONS 

Fig. 1. The total DOS and the partial Mn-d and 0-2p DOS of the monoclinic 
LiMn02 in antiferromagnetic state calculated at optimized lattice parameters. Upper 
most panel shows total DOS. The middle and lower panels show partial Mn-d and 0-2p 
DOS respectively. 

Fig. 2. Partial d-DOS of AFM monoclinic LiNi .25Mn . 75 O 2 . Panel (a) shows d- 
DOS of Ni 2+ . Panel (b), (c) and (d) show d-DOS for Mn ions that are neighbors of Ni 
dopant. One Mn adopts 4+ oxidation state, while other two remain in the 3+ oxidation 
state. 

Fig. 3. The function H(x) at x — 0.25 for dopants Mg, Ni, Co and Fe, calculated 
within GGA+U approximation. 
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